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Abstract. When a finite volume of an etching solution conies in contact with a dis- 
ordered solid, a complex dynamics of the solid-solution interface develops. Since only 
the weak parts are corroded, the solid surface hardens progressively. If the etchant 
is consumed in the chemical reaction, the corrosion dynamics slows down and stops 
spontaneously leaving a fractal solid surface, which reveals the latent percolation crit- 
icality hidden in any random system. Here we introduce and study, both analytically 
and numerically, a simple model for this phenomenon. In this way we obtain a detailed 
description of the process in terms of percolation theory. In particular we explain the 
mechanism of hardening of the surface and connect it to Gradient Percolation. 



INTRODUCTION 

Corrosion of solids is an everyday phenomenon with evident practical conse- 
quences and applications [1]. Recent developments of theoretical tools for the 
study of disordered systems and fractals in the context of statistical mechanics, has 
triggered and outburst of activity in this subject. 

When an etching solution is put in contact with a disordered etchable solid, the 
solution corrodes the weak parts of the solid surface while the hard parts remain 
un-corroded. During this process new regions of the solid (both hard and weak) 
come into contact with the etching solution. If the volume of the solution is finite 
and the etchant is consumed in the chemical reaction, then the etching power of 
the solution diminishes progressively and the corrosion rate slows down. When the 
solution is too weak to coorrode any part of the hardened solid surface, the dynamics 
spontaneously stops. Note that the etchant concentration at the arrest time is not 
zero, and, as shown below, its final value is strictly related to the percolation 
critical point. One of the most interesting aspects of this phenomenon is that the 



final connected solid surface has a fractal geometry up to a certain characteristic 
scale a given by its thickness. This is precisely the qualitative phenomenology 
that has been recently observed in a nice experiment on pit corrosion of aluminum 
thin films [2]. As we show below this fractality falls in the universality class of 
percolation theory [3]. 

THE MODEL 

Recently, a simple dynamical model of etching, capturing the aforementioned 
phenomenology, has been proposed and studied [4,5]. The model is reviewed in 
Fig. 1 and described as follows: 

(1) The solid is represented by a lattice of sites exhibiting random "resistances to 
corrosion" Tj G [0, 1] uniformly distributed. It has a width L and a given fixed 
depth Y (in the case infinite). The etching solution has a finite constant volume V 
(V ^> L) and contains an initial number N et (Q) of etchant molecules. Therefore the 
initial etchant concentration is given by C(0) = N et (0)/V. At any time-step t the 
"etching power" p(t) of the solution is proportional to the etchant concentration 
C(t). Hereafter we take the proportionality constant equal to one without loss of 
generality, and we choose p c < p(0) < 1, where p c is the percolation threshold of 
the lattice. 

(2) The solution is initially in contact with the solid through the bottom boundary 
y = 0. At each time-step t, all surface sites with < p(t) are dissolved, and a 
particle of etchant is consumed for each corroded site. Hence, the concentration of 
the solution decreases with t. Moreover, depending on the connectivity criterion 
chosen for the lattice, m(t) new solid sites, previously in the bulk, come into contact 
with the solution. Note that at the next time-step, only these last sites can be 
corroded by the solution, as the other surface sites have already resisted to etching 
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FIGURE 1. Sketch of the etching dynamics in a square lattice: the sites 2, 3, 5 are etched at the 
first time-step as their resistances are lower than p(0). At the same time the number of etchant 
particles in the solution decreases by 3 units, and a new part of the solid is uncovered. 



and p(t) decreases with t. Calling n(t) the number of dissolved solid sites at time 
t, one can write 

P(t + l)=p(t)-^. (1) 

The corrosion dynamics spontaneously stops at the time-step tf, when all the sur- 
face sites have a resistance r > p(tf). 

The main characteristics of the model are: 

(i) The final value Pf — p(tf) is slightly smaller than the percolation threshold p c . 
The difference \pf — p c \ ^ as (L/V) ap with a p ~ 0.45, when the limits L — > +oo 
and V/L — > +oo are taken in the right way [5]. This fact can be explained in 
terms of percolation theory. Indeed, this theory implies that, for p(t) < p c , the 
probability to have a connected path of solid sites each one with r > p(t), then 
stopping corrosion, is equal to one in the right thermodynamic limit. 

(ii) The global temporal behavior of p(t) can be divided (see Fig. 2-a) into two 
regimes: a smooth regime when pit) > p c , and a critical regime when p(t) ~ p c . 
The duration of the smooth regime is measured approximately by t c , defined by 
p{tc) — Pc- The quantity t c is found to scale with the ratio L/V in the following 
way t c ~ V/L. The duration of the critical regime is given by (tf — t c ), and scales 




FIGURE 2. a- Time decay of p(t) (continuous line), compared to the analytical approximation 
p(t) = po exp{—tL/V) in the smooth regime (dashed line). It refers to simulations with p(0) = 1 
and V/L = 1000 in a triangular lattice. The two lines are in good agreement for p(t) > p c . 
b- Three snapshots of the corroded solid (dark grey) by the solution (clear) at three different 
time-steps: initial, intermediate and final. Black regions represent detached islands. The final 
corrosion front shows fractality up to a characteristic length given by its thickness. 



as (tf — t c ) ~ (V/L) at f with a tf ~ 0.55 (apart from a further linear dependence on 
log L of the scaling coefficient due to the extremal nature of tf [6] ) . 

(iii) The global surface of the solid in contact with the solution displays a pecu- 
liar roughening dynamics (three snapshots at almost initial, intermediate, and final 
time are reported in Fig. 2-b). First of all, during the process, the solution can 
surround and detach finite solid islands. Consequently, the global solid surface is 
composed by the sum of the perimeter of the "infinite" solid, here called corrosion 
front, and the surfaces of the finite islands. In the first smooth regime, the corrosion 
is well directed and the corrosion front becomes progressively rougher and rougher, 
while finite detached islands are quite small. In the second critical regime, the 
correlations revealed by the hardening process become important: the dynamics 
becomes locally isotropic generating a fractal front. Detached islands in this last 
period are quite big with a fractal perimeter too. 

(iv) The final corrosion front is fractal with a fractal dimension Df = 1.753 ± 0.05 
(see Fig. 3- a) up to a characteristic scale a which is the final thickness of the front. 
Beyond this length the corrosion front acquires a one-dimensional shape. Also a 
is found to scale with the ratio L/V in the following way: a ~ (L/V)~ a,T with 
a a ~ 0.57 ~ 1/Df. As shown below, all these properties can be explained in the 
framework of percolation theory and in particular of Gradient Percolation [7]. 

(v) Another important feature captured by this simple model is the progressive 
hardening of the solid surface, due to the corrosion only of sites with a resis- 
tance weaker than the etching power of the solution, which, in its turn, decreases 
with time. This hardening effect is described quantitatively by the behavior of 
the normalized resistance histogram 0(r, t) of the global surface sites. Obviously, 
0(r, 0) = 1 with r G [0, 1], while 4>(r,tf) is given in Fig. 3-b. As one can see, the 
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FIGURE 3. a-Box counting measure of the fractal dimension of the corrosion front for length 
smaller than a. b-Final normalized histogram of the global solid surface: simulations results 
(empty points) are compared with the analytical evaluation (continuous line). 
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final global surface is much more resistant to corrosion than the initial one. 



ANALYTICAL RESULT AND PERCOLATION 

In order to have a mathematical insight in the previous numerical results, and a 
more clear explanation in terms of percolation, it is useful to study the temporal 
evolution of the above defined resistance histogram. Let be h(r,t)dr the number 
of surface sites with resistance in the interval [r, r + dr] at time t ((f>(r, t) is the 
normalized version of h(r,t)). Obviously, the number of corroded sites at time t is 
n(t) = Jq^ drh(r,t). Let us define uo(t) = m(t)/n(t). Then, using Eq. (1), we can 
write the two coupled equations: 

rp(t) 

h(r, t + 1) = h(r, t) - h(r, t)6\p(t) - r] + u(t) / dr' h(r', t) 



o 

, x Io {t) dr'h(r',t) 
p(t + 1) = p(t) - * , (2) 

where 6{x) is the usual step Heavyside function. In the first relation of Eq. (2), 
the second term on the right side represents the number of just corroded sites 
with resistance between r and r + dr, while the third one gives the contribution 
to the new histogram coming from the m(t) new surface sites. All the geometrical 
features of the system related to the behavior of p(t) and h(r, t) are described by 
the coefficient u(t). In order to close this system of equations, we introduce two 
different approximations of u>(t) for the smooth and the critical regimes. One can 
show [5] that, since in the smooth regime n(t),m(t) 1 and fluctuations are small, 
the right approximation is u(t) = \/p(t). Using this relation in Eq. (2), one gets 

p(t)=p(0)exp(-t/r), (3) 

with r = — 1/ ln(l — L/V) ~ V/L. The good accuracy of this approximation in the 
smooth regime is shown by Fig. 2-a. Eq. (3) gives t c ~ V/L in good agreement with 
simulations. In the critical regime, the only simple approximations is a mean-field- 
like relation w(t) = l/p c [5], forgetting fluctuations which are important in this 



regime. Putting this relation inside Eq. (2) one obtains, (p c — Pf) ~ yL/V and 



(tf — t c ) ~ yV/L which are not very far from the simulation behaviors. Finally, the 
histogram h(r, t) obtained by Eq. (2) within these approximations is in excellent 
agreement with simulations (see Fig. 3-b). 

The link between this dynamical etching model and percolation criticality can be 
shown by analyzing more deeply Eq. (3) valid in the smooth regime. From Eq. (3), 
one has t c ~ V/L\n\p(0)/p c ]. This implies that the condition for smooth regime is 
t <V/L. Then in this temporal range one can linearize Eq. (3): 



pit) ~p(0)(l -tL/V) 



Since in this regime the solution advances mainly layer by layer into the solid, one 
can say that the layer etched at time t is at a depth y = t. We can map then 
our dynamical model in a static model defined by the following rule: a general 
solid site belonging to the layer y is replaced by the solution with a y-dependent 
probability p(y) ~ £>(0)(1 — y L/V). This static model, if extended to any value of 
p(y) > 0, is known under the name of Gradient Percolation (GP) [7]. In our case 
the gradient Vp{y) of replacement (etching) probability satisfies Vp ~ L/V . Why 
this mapping can be extended also out of the smooth regime is more clear from field 
theory arguments [8] . The GP model is characterized by the following features: 
the boundary (the corrosion front in our case) between the connected etched region 
and the connected non-etched one has fractal features in the universality class of 
critical percolation. In particular it has a fractal dimension Df = 7/4 coinciding 
with the hull fractal dimension of infinite percolation cluster at criticality [9]. The 
fractality is not extended to the whole corrosion front, but it is restricted to regions 
of size a ~ Vp~ 1 / Df , which, considering Vp ~ L/V, is in good agreement with 
our numerical results. Furthermore, the range of p's on this boundary is given 
by Ap ~ aVp ~ (L/V) 1 ~ 1 / Df ', which is well confirmed by simulation considering 
the quantity p c — pf. Note that from previous scaling relations, one can write 
a ~ Ap~ u , where v = l/(-D/ — 1) = 4/3 is the usual correlation length exponent in 
percolation. This shows that in the dynamical etching model a can be interpreted 
as the correlation length of the corrosion front: that is, the fact that pf is not 
exactly equal to p c introduce a finite correlation length in the problem. 
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